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We consider a spherically symmetric self-gravitating massless scalar field enclosed inside a timelike 
orldtube IxS 3 with the perfectly reflecting w 
small generic initial data evolve into a black hole. 



worldtube 1 x S 3 with the perfectly reflecting wall. Numerical evidence is given that arbitrarily 
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Bizori and Rostworowski [l| recently gave evidence that 
any generic small perturbation of anti-de Sitter (AdS) 
spacetime evolves into a black hole after a time (D(e~ 2 ), 
where e measures the size of the perturbation. They ar- 
gued that this instability is a combined result of the fol- 
lowing two effects. First, the timelike boundary at spatial 
infinity of AdS acts as a mirror at which the waves prop- 
agating outwards bounce back and return to the bulk. 
Second, nonlinear resonant interactions of reflected waves 
tend to shift their energy into increasingly high frequen- 
cies, or equivalently increasingly fine spatial scales, which 
eventually leads to the formation of a horizon. There 
arises a natural question whether this kind of turbulent 
dynamics is specific to asymptotically AdS spacetimes, or 
it can also occur in other situations when dissipation of 
energy by dispersion is absent. Motivated by this ques- 
tion, in this paper we consider small perturbations of a 
portion of Minkowski spacetimes enclosed inside a time- 
like worldtube RxS 3 . Admittedly, this problem is some- 
what artificial geometrically, yet we think that it sheds 
some new light on the results of 

We restrict ourselves to spherical symmetry and choose 
a minimally coupled massless scalar field <p(t, r) as mat- 
ter source. It is convenient to parametrize the general 
spherically symmetric metric as follows 



ds 2 



A_ 



dt 2 + A- l dr 2 + rW . 



(1) 



where r is the areal radial coordinate, dfl is the round 
metric on the unit two-sphere, and the functions A and 
N depend on t and r. We introduce auxiliary variables 
$ := 4>' and II := A~ l N<j> (hereafter primes and dots de- 
note d r and dt , respectively) and write the wave equation 
V Q V Q <i> = in the first order form 
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Einstein's equations G p „ = 8ttG (d^d^cj) — ^9nv(d(f)) 2 ) 
take a concise form when expressed in terms of the mass 
function m 



hr(l — A) (in units where AttG = 1) 
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To ensure regularity at r — we require that 

$(i,r) = 0(r), II(f, r) = O (r 2 ) 
m(t,r) = O (r 3 ) , N(t, r) = 1 + O (r 2 ) 



(6) 
(7) 



where we set N(t, 0) = 1 so that t is the proper time at 
the center. 

In the asymptotically flat situation the above system 
has been extensively studied in the past leading to im- 
portant insights about the dynamics of gravitational col- 
lapse. In particular, Christodoulou showed that for small 
initial data the fields disperse to infinity 2], while for 
large initial data black holes are formed [3j. The bor- 
derline between these two generic outcomes of evolution 
was explored numerically by Choptuik leading to the dis- 
covery of critical phenomena at the threshold of black 
hole formation [4|]. The aim of this note is to see how 
these findings are affected by placing a reflecting mirror 
at some finite radius r = R. In other words, instead 
of asymptotically flat boundary conditions, we consider 
the interior Dirichlet problem inside a ball of radius R 
with the boundary condition R) = 0. The Dirichlet 
boundary condition and the requirement of smoothness 
imply that the coefficients of the power series expansions 
at r = R (here p = 1 — r/R) 



(8) 



fc=i 



N(t,r) = J2 N k(t)p\ m(t,r) =J2 m k(t)f 



k=0 



k=0 



are determined recursively by </>i(i), No(t), and mo(t). 
It follows immediately from equation (|5j) that ma(t) = 
M = const, where M is the total energy. From @ the 
total energy can be expressed as the volume integral 



M 



1 r R 

- / A(§ 2 +n 2 )r 2 dr. 
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Taken at t = 0, the expansions © express the compati- 
bility conditions between initial and boundary values. 

We solved the initial-boundary value problem for the 
system ©-(O numerically using the method of lines 
with a fourth-order spatial finite-difference discretization 
scheme. Time integration of evolution equations was per- 
formed with use of an adaptive, explicit Runge-Kutta- 
Dormand-Prince algorithm of order 5(4). The metric 
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to zero. However, a wave packet which is marginally too 
weak to form a horizon on the first implosion, does so 
on the second implosion after being reflected back by the 
mirror. As in the AdS case, this leads to a sequence 
of critical amplitudes e n for which the solutions, after 
making n bounces, asymptote Choptuik's critical solu- 
tion (see Fig. [1}. To track the steepening of the wave 
packet for very small amplitudes, we follow [l] and mon- 
itor the Ricci scalar at the center R(t,0) — ~2H(t, 0) 2 . 
This function oscillates with approximate period 2. Ini- 
tially, the amplitude stays almost constant but after some 
time it begins to grow exponentially and eventually a 
horizon forms (see Fig. As shown in Fig. the 

time of onset of exponential growth T scales with the 
amplitude of initial data as T ~ £~ 2 , which indicates 
that arbitrarily small perturbations (for which it is im- 
possible numerically to track the formation of a horizon) 
eventually start growing. 



FIG. 1: Apparent horizon radius tah (top) and corresponding 
formation time tAH (bottom) as a function of the amplitude of 
initial data (1121) . At critical points lim + tah{£) = 0, while 
the horizon formation time exhibits jumps of size tAH(£n+i) — 
tAH(£n) ~ 2 (time in which the pulse traverses the cavity back 
and forth). 

functions were updated by solving the slicing condition 
© and the Hamiltonian constraint (Q}. Preservation of 
the momentum constraint ([5]) and conservation of the 
total mass ([9]) were monitored to check the accuracy of 
the numerical scheme. Let us point out that this fully 
constrained evolution scheme is very efficient computa- 
tionally and easy to make parallel because the update of 
metric functions reduces to simple integrations. Namely, 
from ([3]) we have 

log N(t,r) = - I s(<S>(t,s) 2 + n{t,s) 2 )ds, (10) 
Jo 

and from the combination of (j3|) and (jl} we get 

r J a N(t,s) 

Numerical results presented below were generated from 
Gaussian-type initial data of the form (without loss of 
generality we set R = 1) 

$(0,r)=0, n(Q,r) = eexp (-32 tan 2 -rj . (12) 

These initial data vanish exponentially as r — > 1 so com- 
patibility conditions are not an issue. The results are 
very similar to those of [lj , as can be seen by comparing 
Figs. [Hand [2] with the analogues figures in [l]. For large 
amplitudes the evolution is not affected by the mirror; 
the wave packet rapidly collapses, forming an apparent 
horizon at a point where the metric function A(t, r) goes 
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FIG. 2: (a) The function II 2 (t, 0) for solutions with initial 
data (|12p for several moderately small amplitudes. For clar- 
ity of the plot only the envelopes of rapid oscillations are 
depicted, (b) The curves from the plot (a) scaled according 
to e _2 n 2 (et,0). 

In [l[ the numerical results were corroborated by a non- 
linear perturbation analysis which demonstrated that the 
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instability of AdS is caused by the resonant transfer of 
energy from low to high frequencies. For the problem 
at hand, the evolution of linearized perturbations is gov- 
erned by the free radial wave equation cj> — r^ 2 (r 2 (f>')' , 
hence the eigenfrequencies and eigenmodes are 
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As in AdS, the spectrum is fully resonant (that is, the 
frequencies uu are equidistant), so the entire perturbation 
analysis of [l| can be formally repeated in our case. We 
say 'formally' because, in contrast to the AdS case, the 
eigenmodes (| 13[) violate the compatibility conditions at 
r = R and therefore they cannot be taken as smooth 
initial data. 

The transfer of energy to higher modes (which is equiv- 
alent to the concentration of energy on smaller scales) can 
be quantified by monitoring the energy contained in the 
linear modes Ej = Il| + u)j 2 <5&, where = (A 1 / 2 ^, e'j) 
and Ilj = (A 1 / 2 !!, ej), with the inner product defined as 
(/>&) : ~ Jq 1 f{r)g{r)r 2 dr. Then, the total energy can be 
expressed as the Parseval sum M = Y^jLi The 
evidence for the energy transfer is shown in Fig. [3] which 
depicts a Sobolev-type weighted energy norm E(t) = 
J^jLi^Ejit). The growth of E(t) in time means that 
the distribution of energy shifts from low to high frequen- 
cies. The characteristic staircase shape of E(t) indicates 
that the energy transfer occurs mainly during the sub- 
sequent implosions through the center. This observation 
leads to the conclusion that the only role of the mirror is 
to reflect the pulse so that it can be focused during the 
next implosion. 
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FIG. 3: Plot of the weighted energy norm E(t) for the solution 
with initial amplitude e = 8. The steep bursts of growth occur 
when the pulse implodes through the center. 

Another aspect of the turbulent cascade is shown in 
Fig. S] which depicts the spectrum of energy (that is, the 
distribution of the total energy over the linear modes) for 
the solution with initial data (|T2"j) with e = 6. Initially, 



the energy is concentrated in low modes; the exponential 
cutoff of the spectrum expresses the smoothness of initial 
data. During the evolution the range of excited modes 
increases and the spectrum becomes broader. Just before 
horizon formation the spectrum exhibits the power-law 
scaling Ej ~ j~ a with exponent a 1.22 . We note that 
the power-law spectrum with a similar exponent was also 
observed in the AdS case As pointed out in [l|, the 
black hole formation provides a cut-off for the turbu- 
lent energy cascade for solutions of Einstein's equations 
(in analogy to viscosity in the case of the Navier-Stokes 
equation). It is natural to conjecture that the power-law 
decay is a consequence of the loss of smoothness of the 
solution during collapse, however we have not been able 
to compute a analytically. 
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FIG. 4: Energy spectra at three moments of time (initial, 
intermediate, and just before collapse) for the solution with 
initial amplitude e = 6. The fit of the power law Ej ~ j~ a 
performed on the interval j £ [16, 128] gives the slope a w 
1.22. 

Close parallels between the results presented in this 
work and [H,@ indicate that the turbulent behavior is not 
an exclusive domain of asymptotically AdS spacetimes 
but a typical feature of solutions of Einstein's equations 
with compact (or effectively compact) spatial slices. This 
answers the question about the role of the negative cos- 
mological constant A posed at the end of the only 
role of A is to generate the timelike boundary at spatial 
and null infinity. 

Finally, let us point out that instead of the Dirichlet 
condition at the boundary of the cavity U(t, R) — one 
can alternatively consider the Neumann boundary con- 
dition $(i, R) = (for which the energy flux through 
the boundary vanishes as well). In the Neumann case 
the eigenfrequencies of the linearized problem are only 
asymptotically resonant, nonetheless numerical simula- 
tions show similar turbulent behavior. 
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